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TRANSFORMATIONS  OF  GAUSSIAN  LIGHT  BEAMS 
CAUSED  BY  REFLECTION  IN  FEL  RESONATORS 


I.  INTRODUCTION 

Free  Electron  Lasers  (FELs)  operating  as  oscillators^”^  require  the 

8  9 

trapping  of  light  pulses  between  systems  of  mirrors  (resonators).  ’  These 
pulses  are  repeatedly  amplified  via  synchronous  interaction  with  electron 
pulses  passing  through  the  wiggler.  The  radiation  produced  by  the 
stimulated  emission  is  confined  within  a  narrow  cone  along  the  beam  axis. 
Therefore,  the  vector  potential  can  be  represented  as  a  superposition  of 

1  Icz 

Gaussian  modes.  Those  are  the  free  space  eigenmodes  A  (r)  =  e  A  (r)e 

Dll  Dti  inn 

where  is  the  polarization  vector,  of  the  paraxial  equation, 

-  2ik  If  =  0.  (1) 

Equation  (1)  is  the  kj^  <<  k.  =  w/c  limit  of  the  exact  wave  equation.  The 

simplest  oscillator  configuration  is  that  of  an  open  resonator  with  two 

opposed  identical  mirrors.  The  vacuum  eigenmodes  for  this  arrangement  are 

also  expressed  in  terms  of  the  paraxial  eigenmodes.  Their  detailed 

structure  can  be  described  in  terms  of  either  Gaussian-Hermite  functions  in 

rectangular  coordinates,  or  Laguerre  functions  in  polar  coordinates.  In 

both  representations  all  the  eigenmodes  with  given  wave  number  k  are 

1/2 

characterized  by  two  independent  parameters:  the  waist  w  =  (2b/k)  and 

2  2 

the  curvature  of  the  wave  front  1/R  =  z/(z  +  b  ),  where  z  is  the  distance 

from  the  waist  position  and  b  is  the  Rayleigh  length  (Fig.  1). 

The  electron  beam  is  an  optically  active  medium  that  alters  the 

characteristic  parameters  of  the  radiation  after  each  passage.  During  the 

build-up  period  the  modal  content  and  the  structure  of  the  light  pulses 

inside  the  oscillator  will  change.  A  numerical  method  has  been  developed 

recently  optimizing  the  representation  for  the  amplified  radiation.  In  the 

11  12 

source  dependent  expansion  ’  the  waist  size  and  the  curvature  of  the 
elected  modal  basis  is  tailored  according  to  the  driving  source  term.  That 
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minimizes  the  number  of  modes  required  to  describe  the  light  beam.  In 
general,  the  curvature  and  waist  size  of  these  modes  does  not  match  the 
curvature  and  waist  of  the  vacuum  eigenmodes  for  the  resonator.  Therefore, 
the  transfer  matrix  for  a  given  mirror  must  be  known  for  arbitrary  incoming 
modes.  This  need  stems  from  computational  as  well  as  physical  reasons. 

The  knowledge  of  the  cavity  reflection  matrix  R,  together  with  the  gain 
matrix  G  through  the  wiggler,  is  necessary  in  determining  the  potential  for 
steady  state  operation. 

The  study  of  the  reflection  matrix  must  include  the  effects  of 
deflecting  the  light  beam,  in  addition  to  finite  mirror  size  and  curvature 
mismatches.  During  high  power  operation,  grazing  mirror  incidence  may  be 
necessary  to  avoid  exceeding  the  dielectric  breakdown  limit  for  the 
reflecting  surface.  Also,  in  case  of  a  high  per-pass  gain  with  optical 
guiding,  the  spot  size  for  the  reflected  radiation  could  be  much  larger 
than  the  incoming.  In  two  mirror  resonators,  the  reflected  radiation  could 
then  damage  the  wiggler.  Therefore,  ring  resonators,  including  three  or 
more  mirrors,  must  be  employed  for  the  deflection  and  recirculation  of  the 
light  pulses. 

We  are  interested  in  cases  when  the  reflected  radiation  remains 
focused  along  some  direction  z^  making  an  angle  ♦  with  the  incoming  z^. 

Then  the  reflected  vector  potential  will  also  be  expandable  in  free  space 
eigenmodes  Ap^(r^)  of  the  paraxial  equation  in  the  new  direction.  The 
mirror  surface  generating  focused  reflection  in  the  desired  direction  can 
not  be  arbitrary  but  must  be  appropriately  defined.  The  angle  of 
deflection  ♦  will  enter  the  equation  defining  the  mirror  surface.  The 
other  surface  parameter,  namely  the  curvature  1/R|^>  is  a  free  parameter. 

It  determines  the  curvature  l/R^  for  the  outgoing  modes  given  the  curvature 
1/R^  of  the  incoming  modes.  In  case  of  reflection  by  an  arbitrary  surface, 
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the  scattered  radiation  cannot,  in  general,  be  covered  by  the  paraxial 
modes  that  do  not  form  a  complete  set  in  three  dimensions. 

A  single  incident  mode  will,  in  general,  be  partially 

reflected  into  different  modes  Ap^(r^)  where  (m,n)^(p,q).  This  is  caused 
by  the  deflection  of  the  light  beam,  the  finite  size  of  the  mirror  and  the 
curvature  mismatches.  Reflection  into  other  modes  will  affect  the 
interaction  between  the  electron  beam  and  the  radiation  in  a  number  of 
ways.  First,  the  rms  radius  of  the  light  beam  will  change,  affecting  the 
matching  beam  condition.  Second,  the  light  pulse  will  spread  axially 
because  of  dispersion  among  different  modes,  since  the  phase  velocity 
depends  on  the  modal  number  (m,n).  Third,  different  phase  shifts  among  the 
various  modes  during  reflection  may  render  these  modes  out  of  phase  after  a 
number  of  bouncings  off  the  resonator.  For  the  above  reasons  the  fraction 
of  radiation  scattered  into  other  modes  will  contribute  to  the  losses  in 
FEL  oscillators. 

The  method  for  obtaining  the  reflection  matrix  is  outlined  in  Sec.  II. 
The  definition  of  the  appropriate  mirror  surface  is  given  in  Sec.  III.  In 
Sec.  IV  the  integral  expressions  for  the  matrix  elements  are  derived.  An 
analytic  expansion  in  powers  of  a  small  parameter  (of  the  order  of  the 
diffraction  angle)  is  given  in  the  same  section.  Some  limiting  cases  are 
examined  in  Sec.  V.  In  Sec.  VI  the  reflection  of  the  fundamantal  mode 
(0,0)  is  studied  in  detail.  Section  VII  deals  with  cross-coupling  effects 
among  the  vector  components  of  the  radiation. 
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II.  OUTLINE  OF  THE  METHOD 


The  free  space  eigenmodes  of  the  paraxial  wave  equation  have 

the  general  form 


A  (r) 
mn' 


u  (r;W) 
mn^  ’  ’ 

/  2,1/2 


Itz  +  k(x^  y^) 

^  nT\  y  —  V 


2R(2) 


i  S  (z) 
mn'  ' 


(2) 


The  first  exponential  in  (2)  contains  the  rapidly  varying  phase  on  the 

wavelength  scale  X  =  2n/k.  The  wavefronts  are  spherical  with  radius  of 

2  2 

curvature  R(z)  given  by  1/R(2)  =  z/(z  +  b  ).  The  spot  size  of  the 

2  2  1/2  1/2 

radiation  envelope  is  W(z)  =  w(l  +  z  /b  )  ,  where  w  =  (2b/k.)  is  the 

waist,  and  the  distance  z  is  measured  from  the  position  of  the  waist.  The 

amplitude  squared  of  the  mode  drops  by  1/2  over  a  distance  equal  to  the 

Rayleigh  length  b  (also  known  as  confocal  parameter).  Most  of  the 

radiation  is  confined  within  a  cone  parametrized  by  the  diffraction  angle 
1/2 

=  W/z  =  (X/bn)  .  The  structure  of  the  amplitude  profile 

depends  on  the  elected  coordinate  system.  U|^^(r;W)  contains  the  slow 

spatial  variation  equivalent  to  a  small  wave  number  perpendicular  to  the  z- 

direction.  Higher  modes  correspond  to  an  increasing  effective  kj^, 

producing  the  slow  phase  shift  expressed  by  the  term  exp  [iS^^(z)].  For  a 

given  k,  the  mode  is  completely  defined  by  the  two  independent  parameters  R 

and  w  (or  any  combination  of  two  out  of  the  four  quantities  R,  w,  z  and  b). 

The  geometry  of  the  reflection  is  illustrated  in  Fig.  2.  The 

subscripts  i  and  o  denote  the  coordinate  system  used  for  expressing 

incoming  and  outgoing  modes,  r^  is  defined  with  the  z^  axis  along  the 

direction  of  incidence  and  r  has  the  z  axis  along  the  direction  of 

o  o 

reflection.  The  origins  are  displaced  from  the  mirror  center  by  1^.  and  1^ 
respectively,  where  1^  is  the  distance  of  the  minimum  waist  w^,  for  the 


4 


incoming  radiation  and  is  the  distance  of  the  minimum  waist  for  the 

reflected  mode.  A  third  coordinate  system  with  the  origin  at  the  mirror 

center  and  z  aligned  with  the  mirror  axis  will  be  useful  in  the 
m 

computations.  Underlined  quantities  r.,  and  r^  stand  for  the  mirror 
surface  coordinates  in  each  reference  frame.  The  transformations  among  the 
various  frames  are  defined  by 


6 

X.  =  X  cos  f  - 
1  s  2 


sin 


X  =  X  cos  $  + 
os  2 


4 

2’ 


^i  =  ^s’ 


(3a) 


^o  =  ^s- 


(3b) 


z . 
1 


4 

z  cos  ?  +  X 
s  2  s 


sin  I  >  1., 


z 

o 


4 

z  cos  ?  -  X  sin 
s  2  s 


1  . 

o 


We  consider  incoming  radiation  of  given  curvature  and  of  arbitrary 
amplitude  profile  A^(rj),  consisting  of  various  modes  (m,n)  with  the  same 
Rj(z).  If  both  incident  and  reflected  radiation  are  expanded  into 
eigenmodes, 


A^r.)  =  i;  c^ 

m,n 


u  (r . ) 

mn  i' 


mn 


.1/2 


i  S 


mn 


(4a) 


A°(ro) 


=  e 


it  (r  ) 
o'  o' 


u  (r  ) 
,o  pq'  o'^ 


p.q 


pq 


r  z 


iS 


pq 


(4b) 


where 


4i(r) 


(4c) 
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the  relation  among  the  Incident  and  reflected  expansion  coefficients  c^ 

c°  is  formulated  by 
pq 

c°=  R  c\  (5a) 


or 


pq 


=  E  r' 


m,n 


mn 

pq 


mn 


(5b) 


where  are  the  elements  of  the  reflection  matrix  R. 

pq 

We  examine  the  case  when  the  mirror  dimensions  p  are  much  larger  than 

the  wavelength  X,  X  <<  p  (otherwise  diffraction  rather  than  reflection 

would  prevail).  We  also  assume  that  the  angle  C  subtended  by  the  mirror 

C  =  where  parametrizes  the  radius  of  curvature,  is  small,  of  the 

order  of  the  diffraction  angle  0^,  C  ~  6^  ~  e.  The  v-th  component  of  the 

reflected  vector  potential  at  distance  Ir  -  r  I  >>  X  from  the  mirror 
^  '0-0' 

surface  S  is  then  given  by 


A  /  .(r  )  =  - 
(V)  o'' 


ik 

2n 


A®  (r  )  (n  •  6r). 
(V) 


(6) 


In  Eq.  (6)  n-6r  is  the  obliqueness  factor  where  6r  =  (r^-r^)/ |r^-r^ |  and  n 

is  the  normal  unit  vector  to  the  reflecting  surface.  The  surface  element 

ds  is  given  by  ds  =  6[z  -f(x  ,y  )ldx  dy  dz  where  z  =  f(x  ,y  )  is  the 
*  'o  o  o  o  o  •'o' 

surface  equation.  Equation  (6)  is  the  convolution  of  a  source  term  A  (r^) 
at  the  mirror  surface  with  the  propagator  exp(ik |r^-r^  | )/ |r^-r^  | ,  i.e.,  a 
superposition  of  spherical  waves  originating  at  S.  The  source  term  A^(r^) 
is  specified  from  the  incoming  vector  potential  A^(r.)  through  the  boundary 
conditions  and  the  coordinate  transformations  (3).  We  will  assume  a 
perfectly  conducting  surface,  where  the  incident  and  reflected  fields  are 
related  by 
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(7a) 


A®  =  -A^  +  2  (n  •  A^)  n, 

and  n  is  the  normal  unit  vector  to  the  reflecting  surface.  The  second  term 
in  (7a)  introduces  a  coupling  among  different  vector  components,  caused  by 
the  mirror  curvature.  This  cross-coupling  is  small  and  disappears  in  the 
plane  mirror  limit, 


=  - 


(7b) 


(V)  “  "  (N))’ 

i  s 

where  A  and  A  are  expressed  in  the  incoming  and  outgoing  coordinate 
systems  respectively.  Because  of  the  linear  superposition  principle,  Eq. 
(6),  the  cross-coupling  contribution  can  be  added  separately,  and  will  be 
deferred  until  Sec.  VII.  In  the  next  three  sections  we  will  treat  the 
reflected  vector  components  as  independent  scalars,  according  to  (7b),  that 
corresponds  to  a  phase  shift  by  it  during  reflection.  Most  of  the 
computations  will  be  performed  on  the  mirror  surface.  To  simplify  the 
notation  from  now  on,  we  drop  the  bar  (_)  under  the  mirror  coordinates  r. 
Subscripted  quantities  such  as  r^,  r^,  r^  will  signify  the  mirror 
coordinates  in  each  reference  frame.  Unsubscripted  coordinates  will  denote 
the  observation  point  in  the  reflected  radiation  frame  of  reference. 

Ve  seek  cases  when  the  reflected  radiation  propagates  along  z^, 
contained  within  a  cross  section  of  dimensions  x,y  <<  z-z^.  The  expansion 
|r-r^j  =  (z-z^)  {1  +  ((x-x^)^  +  (y-yQ)^|/2(z-z^)^}  replaces  the  full 
propagator  inside  (6)  with  the  paraxial  propagator  in  that  direction, 

A°(r)  =  JJds  A(r^)  (n  •  6r)  U_j^(r,r^),  (8a) 
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where 


(r. 


= 


-ik(z-z  ) 

ik  e _ 

2n  z-z 

o 


-ik 


(x-Xo)^+  (y-y^)^ 


(8b) 


Expression  (8)  is  the  approximation  of  the  exact  solution  (6)  to  order 
2  2  2  2 

((x-x^)  +  (y-y^j)  )/2(z-z^)  ~  c  .  It  is  valid  provided  the  surface  S 

produces  focused  reflection  along  the  desired  direction.  Otherwise  the 
paraxial  limit  will  fail  to  encompass  all  the  radiation  contained  in  the 
original  expression  (6).  The  geometry  of  the  mirror  that  is  compatible 
with  the  above  approximation  will  be  obtained  during  the  computation  of  the 
reflection  matrix. 

It  is  known  that  the  profile  of  a  given  eigenmode  A  (x  ,y  ,z  )  at  z 

^  ®  mn  o''o  o  o 

is  generated  by  the  propagator  U.  (r,r  )  acting  on  the  mode  A  (x,y,0)  at  z 

K.  O  IDD 

=  0.  The  inverse  propagator  U  .  (r,r  )  therefore  reproduces  A  (x,y,0)  from 

“■K  o  nin 

A  (x  ,y  , z  ) , 
mn'  o’-'o  o'* 


u  (x  ,y  ,z  ) 

,  ,  mn  0-^0  o' 

dx  dy^  - V— To  e 

o  o 


1  + 


2ol/2 


2  2 
X  +  y 
,  _o_ 

o^  2R(z  /  ^ 

o'  U_^(r,r^) 


r  X  +  y  , 


u,„(x.y,0).  (9) 


This  suggests  multiplying  and  dividing  the  integrand  inside  (8a)  by 

2  2  1/2 

exp(i#(r^)]  /  (1  +  z^  /b^  J  ,  recasting  (8a)  in  the  form, 


pp  iA(r)  -i$(r) 

^°(r)  =  Jjds  e  °  S(r^)  e  °  °  U_j^(r,r^), 


(10) 


where  the  source  term  S(r^)  is, 


S(r^)  =  A^[r.(r^)]  (n  •  &r) 


2^^  yl/2 

o  '  o' 


(11) 
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and  the  phase  term  given  by, 


A(ro)  =  k 


z .  (r  ) 

i'  0“^ 


+  z 


x.^(r  ) 
1  '  o'^ 


+ 


2 

X  + 

o 


2Ri(ro) 


2R  (r  ) 
o'  o' 


(12) 


The  phase  d(r^)  depends  on  the  angle  <j>  through  the  coordinate 
transformations  Eqs.  (3). 

The  term  exp[iA(r^)J  is  varying  rapidly,  on  the  scale  of  the 
wavelength  X.  Therefore,  its  convolution  with  the  slowly  varying  source 
term  over  an  arbitrary  surface  will  be  vanishingly  small.  In  general,  this 
corresponds  to  radiation  scattering  where  only  a  small  fraction  of  the 
incoming  radiation  is  reflected  along  the  considered  direction  The 
integral  (10)  will  be  finite  only  when  it  is  possible  to  satisfy  the 
condition  ^(r^)  =  constant  over  some  surface  S.  If,  in  addition,  S  is  much 
larger  than  X,  expression  (10)  will  be  finite  only  within  a  narrow  angle  St 
around  t.  This  guarantees  that  the  reflected  radiation  remains  focused 
along  that  direction.  Therefore,  a  condition  that  the  exact  reflected 
radiation  (6)  be  fully  covered  by  the  paraxial  limit  (10)  is  that 

A(r^)  =  constant,  (13) 


along  the  surface  S.  Accordingly,  the  optical  path  is  the  same  along  the 
rays  connecting  an  incoming  wave  front  with  its  mirror  image  (reflected) 
wave  front. 

Requirement  (13)  defines  the  appropriate  mirror  surface  z^=f^(x^,y^; ♦) 
for  reflection  in  the  elected  direction.  Assuming  that  f^  is  found,  we  may 
express  z^  in  terms  of  x^,  y^  and  use  the  constancy  of  Mr^)  over  S, 
reducing  (10)  into 


A°(r) 


a(Xo,yo)  e 


-i^  (x  ,y  ) 

O  0  0 


U_^(r,ro). 


(14) 
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s  expanded  in  terms  of  u^n^^o’^o^ 

form  a  complete  set  in  tvo  dimensions, 


<^(x  ,y  )  =  1]  R™  u  (x  ,y  ;y  ). 
'  o  ■'o'  ^  mn'  o  •'o'  o' 

m,n 


(15) 


The  expansion  coefficients  R  for  Gaussian  incoming  radiation  of 
arbitrary  profile  a(x^,y^)  are  given  by 

r"”^  =  ffdx  dy  a(x  ,y  )  u  (x  ,y  ;W  )/rrdx  dy  u  ^(x  ,y  ;W  ).  (16) 
JJ  o  ■'o  '  o  •'o'  mn'  o’-'o  o'  JJ  o  •'o  mn  '  o’-'o  o'  '  ' 


The  radiation  spot  size  at  the  location  of  the  mirror  center  is  a  free 
parameter,  yet  to  be  specified.  Each  choice  of  generates  an  equivalent 
representation  for  a(x^,y^). 

Upon  substituting  expansion  (15)  inside  the  integral  (lA)  and  using 
the  property  (9)  for  the  inverse  propagator  U  the  reflected  vector 
potential  assumes  the  final  form 


A°(x,y,0)  =  Z  r"""  u^^(x,y;W^),  (17) 

m,n 

?  ?  t/?  1/7 

where  Wq(z)  =%(!■'■  %  “  (2b^/lt)^  .  Expression  (17)  is  a 

complete  decomposition  of  the  reflected  radiation  into  paraxial  eigenmodes 
for  incident  radiation  of  arbitrary  profile.  Therefore,  condition  (13) 
that  defines  the  mirror  surface  is  sufficient  for  the  full  reflection  of 
paraxial  (Gaussian)  incoming  light  beams  into  paraxial  beams  only.  The 
fraction  of  the  electromagnetic  flux  incident  on  the  mirror  is  conserved 
after  reflection.  If,  on  the  other  hand,  (13)  is  seriously  violated,  the 
paraxial  modes  are  inadequate  to  include  all  reflected  radiation,  and  the 
incident  flux  is  not  conserved  by  expressions  similar  to  (17). 
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III.  MIRROR  SURFACE. 


To  obtain  the  equation  for  S  ve  express  all  quantities  inside  (12)  in 

the  mirror  coordinate  frame  applying  the  transformations  (3a)  and  (3b). 

2 

Using  the  scaling  ~  ^  ~  e  we  obtain  from  (13) 


z 


s 


1 

2R  cos 
m 


2 


cos 


2  i 


(18a) 


where 


1  1  1 

R  "  2R  2R.  ■ 

m  o  1 


(18b) 


Equation  (18a)  is  the  analytic  expression  for  a  paraboloid  surface. 

R  parametrizes  the  mirror  curvature,  being  positive  or  negative  for  a 
m 

convex  or  concave  mirror  respectively.  The  surface  is  reflection  symmetric 

with  (zx)  and  (zy)  as  the  symmetry  planes;  there  is  no  rotational 
s  s 

symmetry  around  z^.  Surface  (18a)  can  also  be  approximated,  to  second 
2  2 

order  in  (x  /R  )  ,  (y  /R  )  by  hyperboloids  or  ellipsoids  defined  by 
s  rn  s  ni 

("s  -  1)  -  2  -  2’ 

(==s  ^  1]  ^  2  ^  ^s^  =  2- 

All  the  surfaces  become  spherical  in  the  limit  of  perpendicular  incidence  + 
=  0,  and  plane  mirrors  when  R^  -»  ®.  Using  the  definition  of  the  curvature 
for  the  paraxial  modes,  Eq.  (2),  and  the  fact  that  R  »  b  in  cases  of 
interest,  we  obtain  from  (18b) 


Relation  (20)  defines  the  curvature  of  the  reflected  modes  from  the 
incoming  mode  curvature  and  the  curvature  of  the  mirror. 
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Equations  (18)-(20)  imply  that 


MCg)  a  a[rj(rg),r^(rg)]  =  const.  +  oTkp  1,  (21) 

'^m 

where  p  parametrizes  the  mirror  size.  A  more  complicated  surface  equation 

(higher  than  quadratic  in  x,  y,  z)  is  required  to  improve  the  constancy  to 

a  higher  order.  In  the  next  section  the  reflection  matrix  will  be  computed 

by  expansion  in  powers  of  kp  »  1,  the  approximation 

^(tg)  =  constant  is  satisfactory  for  a  first  order  expansion  as  long  as 

p/R  -  1/kp.  In  case  that  p/R  >  1/kp,  A(x  ,y  )  is  a  slowly  varying 
0)  m  s  s 

function  over  S.  Large  mirrors  require  the  inclusion  of  the  phase  slippage 
term  exp  IiA(x^,yg)l  next  to  the  source  term  a(x^,y^)  in  Eq.  (16). 

The  unit  vector  n  normal  to  the  mirror  surface  is  given  by 


Vf  ♦  ’'s  * 

TWT  =  2  r  ’'s  " 

'  '  m  cos 


cos 

m  / 


z  t 
s’ 


where  f  (x^,  y^,  z^)  is  given  by  Eq.  (18a). 
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IV.  COMPUTATION  OF  THE  REFLECTION  MATRIX 
According  to  the  definition  (5b),  the  R™^  element  of  the  reflection 
matrix  R  is  obtained  from  the  source  term  ®pq(’'o’^o^  inside  (14)  generated 
by  a  single  incident  eigenmode  Apq(rj(r^)].  The  integration  is  performed 
in  the  mirror-aligned  coordinates,  taking  advantage  of  the  existing 
symmetries.  The  coordinates  r^  and  r^,  defining  the  incoming  and  outgoing 
wave  functions,  become  explicit  functions  of  x^,  y^  through  the 
transformations  (3).  The  surface  equation  (14a)  is  used  to  express  in 
terms  of  (x^,yg).  The  mirror  boundary 

Xg  cos  f  +  Yg  =  P  (22) 


is  defined  by  the  intersection  of  the  infinite  surface  (18a)  with  the  plane 
2  2 

z^  =  const  =  2p  cos  (♦/2)/R^.  After  the  above  manipulations,  the 
reflection  matrix  elements  take  the  form 


^mn 

pq 


u  (x  ,y  )u  (x  ,y  ) 

c  ^  c 


s  's 


1 


2  ,1/2 


1  + 


1  + 


-•  ,y  ) 
1  s  s 


b.‘ 

1 


1/2 


iA(X3,ys) 


X  e 


i6^  (x  ,y  ) 
pq'  s’-'s' 


iS  (X  , 

mn'  s 


cos  1  ( 

1  -  ^  sxn 

i 

X  sin 
s  21 

2 

D  2  J 

m 

R 

m 

,(23) 


where 


u  (x  ,y  )  s  u 
mn'  s  •'s  mn 


],  u 


(x  ,y  )  s  u 
pq  s  's  pq 


[xjO 


-ys>’  ^s 


]■ 


(24) 
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2  2 

Expression  (23)  is  correct  to  order  p  /R^  . 

It  will  be  seen  that  R,  as  given  by  (23),  depends  on  four  parameters 

R  H  R(*,  «,  u;  1).  (25) 

♦  is  the  reflection  angle  shown  in  Fig.  2.  a  is  the  ratio  of  the  incoming 

to  the  outgoing  spot  size  at  the  mirror,  a  =  Wj(lj)/W^(l^) .  p  =  p/W^ 

parametrizes  the  mirror  size  compared  to  the  radiation  spot  size.  = 

V  /R  scales  as  the  diffraction  angle  6,  =  W  /I  multiplied  by  the 
om  doo 

curvature  mismatch  R  /R  between  the  mirror  and  the  radiation  wavefronts. 

o  m 

The  spot  size  enters  as  a  free  parameter  because  only  the  curvature  1/Rq 
for  the  reflected  modes  is  specified  by  the  mirror  geometry.  Since  many 
combinations  of  and  1^  apply  to  a  given  curvature  according  to  paragraph 
Eq.  (2),  an  additional  selection  rule  for  is  needed.  Note  that  does 
not  have  to  match  W^.  This  is  obvious  in  cases  when  the  mirror  size  p  is 
smaller  that  W^.  Each  value  of  defines  a  complete  set  of  modes  for  the 
reflected  radiation  and  an  equivalent  representation  for  R. 

Parameters  a,  and  p  can  be  arbitrary.  In  most  cases  of  interest, 
however,  ^  is  small,  t  «  1,  of  the  same  order  as  the  diffraction  angle  9^. 
The  analytic  computation  of  the  matrix  elements  is  carried  out  by  expanding 
the  integral  (23)  in  powers  of  1., 

R  =  R(0)  +  R(l)  +  R(2).  (26) 

Each  representation  of  R  is  tied  to  the  choice  of  the  basis  functions 

u  (r).  The  eigenmodes  u  (r)  are  specified  according  to  the  coordinate 
mn  °  mn 

geometry.  In  the  next  subsections  we  derive  R  in  Gaussian-Hermi te  and 
Gaussian-Laguerre  representations.  For  simplicity,  it  is  assumed  that 
A(x^,y^)  in  Eq.  (23)  is  constant,  i.e.,  kp  ( p/R^)  <<  1. 
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(a).  Gaussian-Hermite  representation 


In  rectangular  coordinates  (x,y,z)  the  functions  U||jj^(x,z;W)  are 


given  by 


2  2 
X  +y 


(27a) 


where  H  ,  H  are  the  Hermite  polynomials  and  a  is  a  normalization  factor, 
m’  n  r  j  » 

setting  the  tot::\l  electromagnetic  flux  carried  by  the  mode  equal  to  unity, 

.-1/2 


/2  /  oHi+n  , 

a  =  -rr  n  2  m!n! 
mn  W  (  J 


(27b) 


The  corresponding  slow  phase  factor  S^^(z)  in  Eq.  (2)  is 


S^^(z)  =  (m  +  n  +  1)  tan~^  ^ - j 


(27c) 


Substituting  inside  (23),  expanding  in  and  performing  the  integrations, 
Eqs.  (23)-(26)  yield 


^mn  ^  ^mn  pq  ^mn 

pq  pq  pq 


(28a) 


^mn 

pq 


«mn  '^pq  ^  ^  Lmn  .  „mnl 

pq  2  I  pq  pqj  ’ 


(28b) 


where  is  a  normalization  factor 


pq 


-1/2 


(28c) 


and  the  phase  is  expressed  by 

*00  "  (b^]  -  ('n+n+l)tan"^  +  k(l.+l  ). 

pq  '•Vj/  Wq/ 


(28d) 
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X  Y 
s  s 


I™  .  JdX  JdY  Bp(aX)  H^(Y)  B„(X)  B„(y) 


-X  -Y 
s  s 


X  Y 
s  s 


C  -  I'**  I'*’'  ® 


S-ii  (X^  * 


(X^  ,  y2) 


,  (29a) 


-X  -Y 
s  s 


(  ,  ,  2  „  ,  ,  H„(aX)  H„(X)  s  ,9  7  9 

I  3  „  1-a  „/„2  „2^  IP  ID  )  /2  ^„Z  „2, 

1^-  2  ^  ^  +  Y  )  - 


(29b) 


-)  ^  (X^  .  y2) 


X  Y 
s  s 


n"'"  =  fdX  [dY—  H  (oX)  H  (Y)  H  (X)  H  (Y) 
PR  J  J  /2  P  9  " 


(X^  +  Y^) 


-X  -Y 
s  s 


.  (29c) 


In  the  rescaled  variables  X  =  cos+/2  /2x  /W  ,  Y  =  /2y  /V  ,  the  surface 

SO  so 

2  2  2  2 

boundary  is  given  by  X^  +  Y^  =  2p  /W^  .  The  lowest  terms  can  be  computed 
directly.  The  matrix  elements  are  computed,  to  first  order  in  in 
Appendix  A  for  large  size  mirror  and  a  =  1. 


(b).  Gaussian-Laguerre  representation 
In  cylindrical  coo 
u^(r,e;W)  take  the  form 


2  2  1/2 

In  cylindrical  coordinates  (r,e,z)  where  tan©  =  x/y,  r  =  (x  +y  )  , 


u«(r,e,B)  . 


'cosp©' 

(dLf  .pf^l 

,sinp9, 

1  W  j  yZj 

1  2r‘ 


(30a) 


where  +p(-p)  signifies  cosine  (sine)  poloidal  dependence,  is  given  by 
1/2 


= 

m 


4  ) 


llW^J 


r.  -1 

l(m+p) !  J 


1/2 


(30b) 
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and  the  are  the  Laguerre  polynomials.  The  corresponding  slow  phase 
S^(z)  in  Eq.  (2)  is, 

8^(2)  =  (2m  +  p  +  1)  tan"^  (30c) 


The  transformations  among  polar  coordinates  representing  the  various 
reference  frames  are 


[2  2  ^  ^s  ^ 

1  -  sin  e^sin  f  "  ^  ®s  2  fj 

111  ^  ^ 


[2  2  ^  ^s  6 

1  -  sin  OgSin  f  ^  ^  R~  2  fJ 


Ola) 


tan 


.  z  sin  f 

cos  f  tao  , 

s  s 


.  z  sin  I 

tan  0  =  cos  ?  tan  0  + - 5-. 

o  2  s  r  cos  0 

s  s 


(31b) 


The  mirror  surface  (18a)  is  expressed  in  polar  coordinates  as 


z 


s 


2  2  ♦  2 
sin  0  cos  i  +  cos  0 
_ s _ 2 _ s 

cos  f 


(31c) 


Applying  similar  computational  procedure  as  in  the  previous  subsection  we 
obtain 


i^^l  A 

^  COS  * 
mn  'mn  mn  2 


rP^  =  cP**  e 


dX  |DP‘>(g^uP‘*(X)+^  jin^  eP'*(X)  [vPJ(X)+i 


(■2a) 
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,1/2 


rPq  -  1-  f— 

mn  2n  Um+p) !  (n+q) !  J  ’ 
and  the  phase  is  expressed  by 

♦nin  =  (2m+p+l)  tan“^  -  (2n+q+l)  tan"^  +  k(lj+l^) 


(32b) 


(32c) 


The  integrals  and  are  given  by 

mn  mn  mn 


2n 

d'^^(X)  =j  de^ 


cos  P  [9^(93)]  cos  q  QqO^) 


1  •  2*  .  2- 

1  -  sin  jp  sin  6 
2  m 


(33a) 


fpq 


2n 

(X)  -I  de. 


sin  0  cos  p 


6.(6^)]  cos  q[e^(e^) 


E  g. 


(1  -  sin^l  sin^e^)^''^ 


UP^(X)  =  (<xV  LP(a^X)  L^(X)  e  ^  , 


mn 


(33b) 


(33c) 


vP^(X)  =  —  -  3  -  X  + 

mn'  "  /2  \  2  2 


[L^U^X)  lJU) 


lP(«X)  lJ(X)J 


2  1 
a  +1 


E  a±i  _ 

(a^X)^  X  ^  LP(a^X)  l‘^(X)  e  ^ 


(33d) 


E  aii 

W  P'’(X)  =  —  (a^X)^  X  ^  LP(a^X) 


mn 


/2 


m<«^>  Ln<">  " 


(33e) 


In  obtaining  (33a)  -  (33e),  X  was  defined  by  X  =  ll-sin^(<}i/2)sin^0)r^/2W^^; 
thus,  according  to  (22)  and  (30),  the  boundary  X^  is  X^  =  2p^/W^^.  The 
lowest  order  terms  for  the  first  few  elements  are  given  in  Appendix  B  for 
arbitrary  deflection  angle  ^  and  a  =  1. 
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V.  LIMITING  CASES 


Vhen  the  mirror  radius  tends  to  infinity  ■^0),  or  in  cases  of 

vertical  incidence  on  the  mirror  =  0),  the  higher  order  corrections  in 
the  reflection  matrix  R  disappear, 

R  =  R(0)  (34) 

in  both  representations.  The  nondiagonal  elements  in  R  stem  from  the 
finite  mirror  size  only.  If,  in  addition,  the  mirror  size  is  very 
large,  u  »  1«  it  is  appropriate  to  take  as  best  representation 

for  the  reflected  radiation.  The  a  =  1  limit  yields 


j^mn  ^mn 

pq  ^  pq* 


(35) 


Thus,  in  case  of  large  curved  mirror  and  vertical  incidence,  or  large 
plane  mirror  and  arbitrary  incidence,  the  reflection  matrix  is  the 
identity  matrix. 

The  case  a  =  1  is  of  special  interest  for  arbitrary  angle  of 
deflection  ♦  and  mirror  curvature  1/R,  as  it  will  be  explained  in  the 
next  section.  For  finite  mirror  size  p  >  W^,  (  u  >  1),  there  exists 
zeroth  order  non-diagonal  terms  inside  R(0).  Since  R(0)  is  independent 
of  the  angle  of  deflection  ♦,  the  finite  mirror  size  yields  the  dominant 
contribution  to  the  reflection  into  modes  different  than  the  incoming. 
The  effects  of  the  deflection  of  the  light  beam  enter  to  first  order  in 
£.,  R(l),  or  higher.  In  the  Hermite  representation  the  elements  Rpq(O) 
couple  mode  combinations  with  m  +  p  =  even,  n  +  q  =  even.  The  elements 
with  either  m  -t-  p  or  n  >  q  odd  vanish  because  of  the  even/odd  symmetry 
of  the  Hermite  functions. 

As  the  mirror  size  becomes  very  large  and  the  limits  of  integration 
in  (23)  are  extended,  the  orthogonality  among  the  various  modes 
becomes  effective.  The  off-diagonal  terms  in  R(0)  become  comparable  to 
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2 

the  first  order  corrections  roughly  when  l/p  ~  ~  e^.  At  the  limit  u 

-*  “  all  the  nondiagonal  elements  of  R  are  reduced  to  order  E,  or  higher, 

Rpq  =  ^  Rpq(l)  +  0(^^),  m  p,  n  ^  q,  (36a) 

and  the  only  matrix  elements  of  zeroth  order  in  f.  are  the  diagonal 


Rmn 

mn 


Rmn 

mn 


(0)  +  0(^^), 


(36b) 


in  both  Hermite  and  Laguerre  represenations.  The  lowest  correction  in  the 

2 

diagonal  elements  is  of  second  order  ^  ,  while  the  first  order  contribution 
disappears.  This  is  consistent  with  flux  conservation  during  reflection  in 
case  of  large  mirror. 

In  obtaining  Eqs.  (28)  and  (32)  it  was  assumed  that  6(x^,y^)  is 
constant  over  S.  According  to  (21)  the  variation  of  L  is  parametrized  by 
J.*  =  (kWj^/R^)  When  (kW^^/R^)  >1,  ^  becomes  of  order  ^  and  the 
effects  of  the  slow  phase  slippage  must  be  retained  in  (23).  This  effect, 
known  as  spherical  aberration,  causes  additional  corrections  R*(l),  of 
order 


R  =  R(0)  +  £.  R(l)  +  R  (t)  +  . 

Spherical  aberration  does  not  disappear  at  the  limit  of  large  mirror 

size,  as  opposed  to  the  effects  discussed  so  far.  In  fact,  when  £,*  >  £., 

it  places  a  lower  limit  on  the  off-diagonal  terms  in  the  reflection 

matrix. 


Rmn  ^ 

pq  -  pq  * 

Perfect  reflection,  requiring  £.  =  0,  is  possible  only  for  plane  mirror  (R^ 

-»  “)  of  large  size. 

The  superposition  principle  can  be  used  to  describe  reflection  from 
more  complex  mirror  surfaces.  In  case  of  a  mirror  with  a  hole  the  surface 
integral  (14)  over  is  expressed  as  Jg  =  -  Jg^  where  is  defined  by 

the  mirror  exterior  boundary  and  S2  is  the  surface  of  the  hole.  The  total 
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reflection  matrix  R  is  given  by  R  =  R(Sj)  -  R(S2),  the  difference  in  the 
reflection  matrices  associated  with  mirrors  and  $2  respectively.  The 
transmission  matrix  T  through  a  screen  with  an  aperture  of  area  S  is  given 
by  T  =  -  R,  R  being  the  reflection  matrix  for  a  mirror  matching  the 
aperture  S.  The  transmission  matrix  for  radiation  diffracted  behind  a 
finite  size  mirror  is  given  by  T'=  1  -  e^*^  R  where  1  is  the  identity 
matrix. 
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VI.  REFLECTION  OF  THE  LOWEST  ORDER  MODE 

The  computation  of  all  the  truncated  Integrals  for  finite  mirror 

surface  is  nontrivial.  Most  applications,  however,  involve  the  (0,0) 

lowest  order  mode  as  the  dominant  mode  in  both  incoming  and  reflected 

radiation.  The  strategy  here  is  to  compute  the  element  R^^  of  the 

reflection  matrix  first.  Then  the  waist  for  the  reflected  modes  W  can  be 

o 

selected  so  that  it  maximizes  Rqq*  The  optimum  representation  condition 


dR, 


00 

00 


3a 


0, 


(37) 


puts  the  maximum  amount  of  the  reflected  radiation  in  the  lowest  order  mode 
(a  different  mode  and  matrix  element  may  be  chosen,  if  desired).  It  is 
pointed  out  that  (37)  does  not  improve  the  properties  of  the  reflected 
radiation.  It  enables  one  to  choose  the  best  representation  in  terms  of 
minimizing  the  coefficients  of  the  undesired  modes  for  the  scattered 
radiation.  Once  is  fixed  by  (37)  then  the  exact  location  and  size  of 
the  waist(s)  for  the  reflected  modes  is  determined  by  solving  the  system  of 
equations 

R“  "  .2.2’ 
o  1  +b 
o  0 


W  =  w 
o  o 


2,  1/2 


1  + 


The  element  R^q  is  identical  in  both  representations 
order  mode  Uqq  is  the  same  in  rectangular  and  cylii.drical 
Performing  the  integration  (29a)  yields  R^q  to  first  order 


(33b) 


since  the  lowest 
coordinates . 
in  f; 
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(39) 


Note  that  the  first  order  term  vanishes  and  the  lowest  correction  is  of 

2 

second  order  in  .  The  exact  dependence  on  the  mirror  size  p  is 

parametrized  by  y  =  while  a  =  parametrizes  the  ratio  of  the 

incoming  and  scattered  radiation  spot  sizes  at  the  mirror.  The  optimization 
condition  3RQQ(0)/3a  =  0  yields,  =  1  +  expl-(l  +  a^)y^] [2M^a^+(2M^+l)a^-l] . 
In  case  that  the  mirror  cross  section  is  much  larger  than  the  spot  size  of 
the  incoming  mode,  y  >>  1,  a  -»  1  and  the  reflected  spot  size  at  the  mirror 
matches  the  incoming, 

Large  mirror  size  is  desired  to  maximize  the  total  reflection 

coefficient.  For  incoming  radiation  of  unity  electromagnetic  flux 
=  |c^|^  =  E  total  reflection  coefficient 

equals  the  reflected  flux  P^, 


I  R 


i  |2 

c  1 


E  E  |r; 

mn  pq  ^ 


(40) 


In  Fig.  3  we  plot  for  the  lowest  order  incoming  mode  as  a  function  of 

y'  =  cos(<fr/2)  =  cos(^/2)  y.  y'  parametrizes  the  size  of  the  mirror 

projection  into  the  plane  perpendicular  to  the  incoming  radiation 

-4 

direction.  The  incoming  radiation  has  a  wavelength  X  =  ly  (10  cm), 
-1  2 

waist  w.  =  2x10  cm  at  distance  1.  =  1.8x10  cm  from  the  mirror  and 

1  1 

3 

radius  of  curvature  (at  the  mirror)  R.  =  8.95x10  cm.  The  mirror  has  a 

3 

radius  of  curvature  R  =  8.95x10  cm,  yielding  reflected  modes  of  R  = 

m  1  j  b  Q 

3 

8.95x10  (again  1^  and  w^  depend  on  the  choice  of  W^).  In  Fig.  4  we 
plot  the  magnitude  of  the  reflection  coefficients  |Ppq|  of  the  lowest 
order  mode  (0,0)  into  the  firstr  25  modes  (p,q)  with  p  <  q  <  5,  as  a 


function  of  y'.  The  deflection  angle  is  90  and  the  ratio  of  the  spot 
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sizes  is  1.  Increasing  mirror  size  maximizes  the  diagonal  element  and 
minimizes  scattering  into  other  modes.  The  spherical  aberration  was 
retained  inside  (23)  in  evaluating  the  matrix  elements.  Its  effect  is 
small,  since  for  the  above  parameters  =  0.28^,  and  a  good  agreement 
is  observed  with  the  constant  A  theoretical  limit.  In  particular,  the 
dominant  off-diagonal  terms  couple  the  (0,0)  incoming  mode  to  the  (1,0), 
(3,0)  and  (3,2)  reflected  modes  only,  according  to  the  selection  rules, 
Eqs.  A(10).  Comparing  Figs.  3  and  4  with  the  next  plots  shows  that  the 
relative  mirror  size  to  the  radiation  spot  size  is  the  most  important 
parameter  to  determine  the  reflection  into  other  than  the  incoming 
modes. 

In  Fig.  5  we  fix  the  mirror  size  p'  =  2  and  the  angle  <(>  =  90°  and 
vary  the  spot  size  ratio  a.  The  best  representation,  maximizing 
and  minimizing  is  obtained  at  a  =  1.  However,  for  small  mirror  u'  = 
0.66,  the  maximum  for  Rq^  occurs  at  a  =  0.70  (see  Fig.  6).  Radiation 
reflected  off  mirrors  smaller  than  the  incoming  spot  size  is  best 
described  by  outgoing  modes  of  reduced  spot  size  <  W^.  Also  note 
from  Fig.  6b  that  for  small  mirror  size  the  total  power  reflected  into 
the  first  25  modes  never  exceeds  80%  of  the  incoming  flux;  even  with 
many  more  modes  remains  less  than  1.  In  Fig.  7  the  reflection 

coefficients  Rp^  are  plotted  as  functions  of  the  angle  of  deflection  <() 
for  fixed  a  =  1,  p'  =  2.  It  is  seen  that,  for  sufficiently  large 
reflecting  surface  and  good  choice  of  the  spot  size  W^,  the  reflection 
matrix  is  not  very  sensitive  to  and  the  off-diagonal  terms  remain 
small . 

The  main  conclusions  so  far  are  summarized  as  follows.  When  the 

mirror  size  is  >  2.5  times  the  incoming  spot  size,  the  fraction  of  the 

2 

incident  power  scatterd  into  different  modes  is  of  order  for 


kW^./R  <  1,  or  (^*)^  for  kW^./R„  >  1.  This  holds  for  a  wide  range  of 
deflection  angles  It  will  be  shown  in  the  next  section  that  cross¬ 
polarization  effects  are  of  the  same  order.  In  most  applications  both  ^ 
*  -2 

and  K  Bre  less  than  10  .  To  this  end,  scattering  losses  will  be 

smaller  than  the  losses  caused  by  the  finite  reflectivity  (i.e., 
absorption)  by  the  mirror,  for  most  dielectrics. 
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VII .  CROSS-POLARIZATION  EFFECTS 

The  curvature  of  the  mirror  surface  produces  a  cross-coupling  between 
the  transverse  components  of  the  incoming  and  reflected  radiation. 
Inserting  expressions  (21)  for  the  normal  unit  vector  to  the  mirror  inside 
the  boundary  conditions  Eq.  (7a),  the  full  source  term 

A®  =  (A  A  A  ®)  for  an  incoming  wave  A^  =  (A  A  0)  is  given  by 
X  y  2  X  y 


A  ®  =  -  A  ^ 


^  +  2  tan^  ^  A 


A  ^  +  2  tan^  5^  A 
y  2  R  X  ’ 


A  ^  =  2  cosf  ^  A  ^  +  2  ^  . 

2  Z  K  X  l\ 

m  m 


1  s 

In  the  above  relations,  the  components  of  A  and  A  are  given  in  coordinate 
systems  aligned  with  the  incoming  and  outgoing  radiation,  respectively. 
According  to  (41)  the  reflection  of  a  plane  polarized  wave  generates 
components  polarized  in  every  direction,  including  A^.  These  cross 
polarization  effects  enter  to  order  E,  and  result  in  a  small  rotation  of  the 
polarization  angle. 

The  radiation  steming  from  the  A^  component  will  propagate 
perpendicularly  to  the  direction  of  interest  z^  and  escapes  the  resonator 
as  pure  reflection  loss.  The  relation  between  the  incoming  and  reflected 
transverse  components,  including  cross-polarization  effects,  is  now  given 


,Q  R)  ry 


The  matrix  3  has  been  computed  in  the  previous  section.  Substitution 
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of  the  additional  cross-terms  in  Eq.  (41)  inside  the  propagator  integral 
(6)  yields 


(xy)pq  (yx)pq  ^pq 


f 


where 


^mn 


2  tani  j - - - - 


1  + 


z  (x  ,y  ) 
o  s  s 


1  + 


"i  ^’'s’ys> 


b.' 

1 


1/2 


X  e 


iS^  (x  ,y  )  -  (x  ,y  ) 
pq  s  •'s  mn  s  •'s 


■ 

X  . 

2x24, 
X  sin 

cos  ^  ^ 

^  -  r  2  - 

s  2) 

D  2  J 

■ 

m 

R 

m 

In  Gaussian-Hermite  representation,  we  obtain 


(43) 


mn 


q'""(1)  =  ’/ll  o'""  e  tan  |  g"*"  , 
pq  pq  2  pq  ’ 


wi  th 


X  Y 
s  s 


GpJ  =  JdX  JdY  Hp(aX)  Hq(Y)  Y  H|^(X)  H^(Y)  e 


2-ii(X^+Y^) 


-X  -y 
s  s 


(44a) 


(44b) 


In  Gaussian-Laguerre  representation,  we  have 
QW(l).^{aiat  j'dXGM(X)  bPOcX), 


(45a) 


where 


1  a+1 

gP^(X)  =  (a^X)^  X  ^  LP(a^X)  L^(X)  e 
mn  m  n 


2  1 
a  +1 


(45b) 
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and 


n  ff  r  n  r  ' 

-  cose  cos  p0.(0)  cos  qe.(0) 

bP^(X)  =  de  - ^ ^ i  .  (A5c) 

2*2 

o  1  -  sin  J  sin  0 

In  both  representations,  cross  polarization  effects  enter  to  order  f;.  In 
case  of  vertical  incidence  (*  =  0)  with  arbitrary  curvature  or  plane 

mirror  (I,  ~  1/R^  =  0)  and  arbitrary  incidence  *,  Q  goes  to  zero. 

Transverse  vector  components  are  reflected  independently  of  each  other  in 
these  two  limits.  Some  of  the  elements  of  Q  (in  both  representations)  are 
given  in  Appendix  C  for  large  (p  »  W.)  mirror. 
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Appendix  A.  Computation  of  the  Hermitian  Matrix  Elements. 

The  integrals  (29)  will  be  evaluated  here  in  case  the  mirror  size  p  is 
much  larger  than  the  incoming  mode  spot  size  Wj,  p  cos  ♦/Z  »  Then  the 

limits  of  the  surface  integrals  can  be  extended  to  inf ini ty^  and  the  spot 
size  for  the  outgoing  modes  matches  that  of  the  incoming  at  the  mirror, 
i.e.,  a  =  1.  Ve  use  the  notation 

=  e  2  H^<X),  <A1) 

the  recurrence  relation 

H^'(X)  =  2n  H^_^(X),  (A2) 

and  the  orthonormality  properties 


00 

—  00 


00 

J  dX  D;^(X)  X  ♦JX) 
—  00 


<X>  =  /n  f2""^n!S  ,  +  2"(n+l)!6  ,1, 

m,n  (  m,n-l  '  m,n+l)’ 

(AA) 


00 

f  dX  (X)  X^  ♦  (X)  a  <X^> 

J  '  ra,n 

—  00 


=  /n  j2"~^n!6  „  +  2"  ^(2n+l)n!S  +  2"(n+2)!6  V,  (A5) 

’  m,n-2  '  m,n  '  '  m,n+2f'  '  ' 


to  obtain 


— fnn  _  ■  i  c  o 

I  =  n  2  n!m!  6  6  « 

pq  P»ni  q,n’ 


(A6) 


m"*"  =  -  I  <X>  f/n  2"njU  , 
pq  2  p,ml.  )  q,n’ 
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-  ^  [2k  2%!  8p_i  „  -  2in  '  Sp,,_l]<^^>q,n’ 

.  i  (P»q*»rr]  <»p,m(''"2""'^)‘q.n- 

pq  ^2  ^  i  i  0  o 

Inserting  expressions  (A3)-(A5)  into  Eqs.  (A6)-(A8)  we  obtain 


'^pq  ^P.m  ^q.n’ 


= 

pq 


(AlO) 


f2n>-l  2n+l')_l/21j 

-  (—  *  ~rr  J*P.m-i 


1  f(;nn^l)(m+2)(m+.: 

/2  ^  « 


p,m+3  ^2 


*  V  8  )  p,n>-3  J  q,n 


Wi-©  ‘p.- 


LlltUillill]  „ 

4  )  q»n+2 


^  [(m+n+l)j^  +  [(V")  ^p,m-l  (  2  )  *p,m+l]^q,n  j‘ 

✓2  *•  i  i  o  o 
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R  is  diagonal  to  zeroth  order.  The  lowest  order  correction  couples  m  with 

m±l,  m+3  in  the  X-direction  and  n  with  n,  n±2  in  the  Y-direction.  The 

reflection  matrix  is  not  symmetric,  R®”  ^  R****.  Also,  it  is  not  invariant 

pq  mn 

' 

to  interchanging  X  and  Y.  This  means  that  the  modes 
with  m  ^  n  are  reflected  differently. 
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Appendix  B.  Computation  of  the  First  Hatrix  Elements  in  Laguerre 
Representation. 

Representation  using  Gaussian-Laguerre  modes  may  be  advantageous  in 

numerical  simulations  because  fewer  Laguerre  modes  than  Hermite  modes  are 

required  to  represent  close-to-axisymmetric  radiation  profiles  with  the 

same  accuracy.  However,  the  computation  of  Eqs.  (33a)  to  (33e)  is  not  so 

straightforward.  The  integrations  (33a)  and  (33b)  for  and  over  the 

polar  angle  0^  involve  trigonometric  functions  of  complicated  arguments 

0j(0g)  and  6Q(0g)*  given  implicitly  by  Eq.  (31b).  There  is  no  simple 

recurrence  formula  for  this  calculation.  The  first  few  elements  are 

computed  here  by  expansions  in  powers  of  <  C  -  Substituting  from 

(31b)  inside  (33)  and  renormalizing  r  ^/2W  ^  =  X/( l-sin^4>/2sin^6) ,  one 

3  O 

obtains,  to  first  order  in  I,, 


1  •  2  i 

1-sin  f 


•  2^ 
sin  9 


(BO) 


D^^(X)  =  -I°^(X)  =  5  X  tan|  f  de  - (Bl) 

2/2  '*  .  2  *  .  2^1 
o  ll-sin  ^  sin  9j 


d“(X)  =  j 


2q 

cos  9 _ 

ll-sin  J  sin  91 


D~^"^(X) 


•  2,, 

sin  0 


^1-sin^  ^  si 


in^e)' 


D~^®(X)  =  -D®  ^(X)  =  -  C  X  tan  I 

2/2  ^ 


2n 

1  '*« 

o 


^1-cos^  ^  tan^9jcos^0 

_  ^ 

(^  •  2  *  •  2  -V 

ll-sin  ^  sin  0j 

(B4) 
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D  ^^(X)  =  D^~\x)  =  0. 


We  only  need  to  zeroth  order  in  obtaining 


E  ^®(X)  =  -  E®~^(X)  =  J 


2a 


de 


4  2 

cos  J  sin  0 


(l-sln^  I  sin^ef' 


(B5) 


(B6) 


and 

eP*!  =  0  +  0(£.)  for  (p,q)  *  (-1,0),  (0-1).  (B7) 

The  integrals  (Bl)  -  (B7)  are  evaluated  using  the  formula 


pn/2 

J 

o 


M  V 

sin  X  cos  X 

(^  1  2  .  2  1  p 

1-k  sin  xj 


(B8) 


where  B(p,q)  =  r(p)r(q)/r(p+q) ,  T  is  the  factorial  function  and  F  is  the 
hypergeometric  function.  The  radial  integrations  for  U,  V  and  W  are 
performed  directly,  using  the  expressions  lP(x)  for  the  Laguerre  functions 
and  the  identities 


r”  -X  -1/2 . 
e  X  dx  =  vn, 

o 

J  e”*  x”dx  =  n!, 
o 

OP 

r  -X  n+l/2,  13  ,  , 

J  e  X  “  2  ’  2  • '  •  ("+l/2)/n.  (B9) 

o 

Again,  we  extend  the  limits  of  integration  to  infinity  assuming  p  cos(^/2) 
»  and  a  =  1.  The  zeroth  order  contribution  is  given  by 

'm'O'  ■ 

rJJ(O)  .  I, 
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R'J  }(0)  =  cos^  I  F(2,|,2,sin2  |), 
Rj|(0)  =  cos  I  F^2,|,2,sin^  |j, 


rJJ(O)  =  cos  I  F(2,|,2,sin^ 

R'J“J(0)  =  cos^  I  F(2,|,2,sin2 

The  first  order  corrections  in  are  given  by 

,1/2  , 


■  -C*''  -  is 


?  Sin  ^  F^2,|,2,  sin^  |j, 


f  FI 


d01/ix  9  .  4 

Rqi^I)  =  °  bJ  16  2 


®  2  Sin^  I], 


Ro?(l)  =  rJ2<1)  =  O' 


1/2 


0'>  ■  -"SI'''  ■  ©  re  I  !)■ 


r}J(1) 


<(1)  - 


I2J  32  2 


2  F^2,^, 2,  sin^  | 


9- 


and 

R-Jg(l) 


Sin  I  t 


.1/2 


rnr  3  „r.  1  „  .2  4')  2  4 

2)  8  2)  -  2 


1  lRg,|,2,  sin^  I) 


±  j  COS  ^  F 


-4.5  +  12'^^^'’ 


L-3.5  -r  i2 


-1/2 


-10 
^  01' 


(1)  -  sin  2  I  -  (f)  fe  I]'  2 


ll 


±  I  cos^  I  Fg,|,2,  sin^  |)[2  -  i 


(BIO) 


(BID 


3  „  ■  2 

,2»2,  sin  D 


2,|,2,  sin^  Ij 
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R  }J(1)  =  sin  I 


i  i  cos^  *  f@.|>2'  2)  [3  5  ^  J2-1/2  }• 


sin^  4)-  cos^  I  f^2,|,2, 


(B12) 


The  (-)  sign  and  the  lowest  row  inside  the  last  square  bracket  in  (B12) 
correspond  to  exchanging  indices. 


mn  nm 

Note  also  that  the  elements  R^~^,  R~^^,  coupling  sine  and  cosine  modes,  are 

ran  mn  ’  ^  °  » 

2 

of  order  E,  or  higher  for  every  m,  n, 


-11  2 


(B13) 
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Appendix  C.  Computation  of  the  Cross-Polarization  Matrix  Elements 

We  compute  here  some  of  the  first  order  cross-polarization  matrix 
elements  in  case  of  large  mirror  size  p  »  and  a  -»  1.  In  the  Gaussian- 
Hermite  representation  we  find  from  (44b),  using  the  notation  of  Appendix 
A, 


g'""  =  2%!  /n  <Y>  & 


q,n  p,m 


=  n  2"'m!  12"  ^n!  S  ,  +  2"  (n+1)!  S  S 

I  q.n-1  q,n+lf  p,m, 


yielding 


q"“‘  (1)  =  t  sin  t  i/n  6  ,  +  /n+1  6  i  ^  S 

pq  2  ]  q,n-l  q,n+lf  p,m. 


(Cl) 


(C2) 


In  Gaussian-Laguerre  representation  we  only  have  to  compute  B^^(X), 
Eq.  (45c),  to  zeroth  order  in  1..  Applying  the  methods  of  Appendix  B,  we 
find 

bP'I  =  0(0,  if  p,q  *  (1,0),  (0,1), 


Noting  that  G^j^(X)  is  the  same  as  Eq.  (33e),  and  inserting  (C3)  and 

(33e)  inside  (45a),  we  obtain 


^00  =  ^00  =  £  sin  I  Pg,  i  2,  sin^  |), 

Ojjd)  =  ojjd)  =  -K  sin  I  f(|,  2,  sin^  |), 

q}J  (1)  =  qJ}  (1)  -  /2  I  sin  I  f(|,  i-,  2,  sin^  (C4) 
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Appendix  D.  Small  Aperture  Limit 

We  have  seen  in  Sec.  V  that  in  case  of  mirror  surface  S^^  with  an 
aperture  of  area  S2  the  reflection  matrix  is  given  by 

R(S;W^)  =  R(Sj;W^)  -  (Dl) 

In  case  that  S2  «  S^^  the  spot  size  W^  optimizing  the  representation 
for  the  scattered  radiation  will  be  determined  predominantly  by  the  surface 
Sj.  Thus,  the  formula  (23)  with  given  from 

3R(S, ;a) 

ai  --  =  0-  (D2) 

can  be  used  for  the  modal  decomposition  of  the  scattered  radiation. 
According  to  Eqs.  (28a)  and  (29a)  for  the  Hermite  representation,  and  Eqs. 
(32a)  and  (33c)  for  the  Laguerre  representation,  the  lowest  order 

2 

contribution  from  a  small  aperture  (>2  scales  as  R(S2;W^)  ~  ^  . 

In  some  cases,  however,  it  is  important  to  know  the  total  radiation 
diffracted  through  a  small  hole,  rather  than  the  modal  decomposition.  In 
case  of  small  apperture  P2, 

P2^«  Igk"^  or  X  »  (D3) 

where  1^  ~  z  is  the  distance  of  the  observation  point  from  the  mirror,  the 
paraxial  approximation,  Eq.  (8b)  is  taken  one  step  further,  setting 

iV  [Ho)  "  Ho)  ]  =  ^  -  2yyo]* 

o  o 
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Substituting  (D4)  inside  (8a)  we  obtain  the  "far  field"  limit  of  the 
diffracted  radiation 


A°(x,y,z=0)  =  ^Jdx  Jdy^A"(x^,y^)(n-8r) 
o 

( -  Ir  -  ^y^o]}’ 


also  known  as  Fraunhoffer  diffraction.  The  condition  (D3)  can  only  be 

valid  for  appertures  much  smaller  than  the  spot  size  at  the  mirror,  p2 

<<  W^,  since  for  P2  ~  (D3)  is  violated,  kp2^  ~  ~  kw^^  (1+1^^/b^^) 

~  b  (1+1  ^/b  ^)  >  1  .  Neglecting  terms  of  order  kx  ^/1  ~  kp-^/1  <<  1 

o  oo  o  ®  mo  Zo 

means  that  terms  of  order  x  /I  ,  x  /I  «  1/kx  ,  where  kx  >1,  must  also 

s  m  s  o  s  s 

be  neglected.  The  source  term  can  be  written  as  A^[x^,y^]  =  A^[x^,yg]. 

Rescaling  variables  to 


,,  <|>  kx 

Kx  =  cos  f  5-, 
o 


K  =  ^ 

y  ^o 


we  obtain 


2  2 

i(pqQ+  k 


A  (x,y,0)  =  e 


°  2  2^  HJ'^ys'^'(’^s,ys> 


iK  X  +iK  y 
g  X  s  y-'s 


According  to  (D7)  the  outgoing  radiation  is  the  Fourier  transform  of  the 

incoming  radiation  in  respect  to  K  ,  K  .  Defining  the  "polar"  coordinates 

X  y 

2  2  1/2  -1  i 

K  =  (K^  +Ky  )  ,  e  =  tan  (K^/Ky),  we  obtain,  for  A 

=  £  u  (x  ,y  ),  the  scattered  radiation 

K,1  pq  pq  s-'s 

2  2 

.  ika  . ,  X  +y  ) 

A°(x,y,0)  =  E  c^  '  2lJ  X 

p,q  o 
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cos 


2  2  i  2 

2  ^  ^s 

W.2 

1 


n 

Y.  Y  j  (Kx  )  J  (Ky  )  i"^"  e  ^ 

^  ^  n  s  m  •'s 

m=-®  n=-“® 


(D8) 


The  zeroth  order  contribution  in  x  /W.  «  1  yields 

SI 

2  .2. 


"00  2nl 


ika..  ik^l_+ 
e 


X  +y 


.O,  i  00  Vo  21, 

A  (x,y,0)  =  c^. 


2  2  *  2 
.  cos  7  +  y 


2  "  ^s 


X  cos  I  JdxJdy^  J^(Kxp 


W.' 

1 


(D9) 


The  waist  size  for  the  Fraunhoffer  inodes  is  given  by  the  zeros  of  the 
Bessel  functions 


kw, 

K(Wf)  ~  ^2  ~ 

o 


(DIO) 


Therefore,  the  diffraction  angle  9^  is 


'^f  X 

B  .  ~  K 

^  \  P2 


(Dll) 


The  requirement  9^  <<  1  for  the  validity  of  the  paraxial  approximation  puts 


a  lower  limit  in  the  aperture  size  P2 


Pj  »  X. 


(D12) 
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In  case  the  aperture  size  is  of  the  order  of  the  vavelength  X  the 
scattered  vavefunctions  are  spherical  rather  than  Gaussian.  Because  the 
overall  effect  of  a  scatterer  with  size  P2  ~  X  is  very  small,  the  familiar 
from  quantum  mechanics  Born  approximation,  involving  perturbation  theory, 
is  applicable  in  that  case. 
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Figure  2  Reflection  geometry. 


44 


Figure  3  Plot  of  the  total  reflection  coefficient  for  the  lowest  order 

mode  as  a  function  of  the  mirror  size  \i'  for  ^  =  90°.  The 

radiation  has  wavelength  X  =  lO'^cm,  waist  w^  =  2xl0"^cm  at 

2 

distance  1.  =  1.8x10  cm  from  the  mirror  and  radius  of  curvature  R 
=  R„  =  8.95xl0^cm. 
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